Colossal Magnetoresistance is a Griffiths Singularity 
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It is now widely accepted that the magnetic transition in doped manganites that show large 
magnetoresistance is a type of percolation effect. This paper demonstrates that the transition should 
be viewed in the context of the Griffiths phase that arises when disorder suppresses a magnetic 
transition. This approach explains unusual aspects of susceptibility and heat capacity data from a 
single crystal of LacrCacsMnOs. 



The term "colossal magnetoresistance" (CMR) has 
commonly been used to describe the very large, magnetic- 
field driven changes in electrical resistivity in oxides 
based on LaMnOa near their second-order, ferromagnetic 
transitions. The largest CMR effects are accompanied by 
other anomalies in magnetic and thermodynamic proper- 
ties. Among these are the failure of the magnetic corre- 
lation length to increase strongly as the transition tem- 
perature Tc is approached from above, the persistence 
of a well defined spin-wave dispersion close to the tran- 
sition [Ql, and an unusual shift in the heat capacitypeak 
to higher temperatures in applied magnetic fields. |^ An 
explanation of this unusual behavior of the heat capacity 
in the context of a Griffiths singularity |^ is the focus of 
the present paper. 

There is now general consensus that the CMR transi- 
tion is a type of percolation in which, due to the double- 
exchange process [Q, bonds become metallic as neigh- 
boring spins tend to align. The strength of the CMR 
effect (along with the transition temperature) depends 
strongly on the ionic size and concentration of the di- 
valent atom that substitutes for La in LaMnOs. The ef- 
fect is nearly absent for La2/3Sri/3Mn03 (Tc = 360 K) 
but quite strong in the present sample Lao.rCao.sMnOa 
{Tc = 218 K) In the low temperature metalHc phase, 
the exchange interactions, as measured by the spin-wave 
stiffness, are the same in Sr- snd Ca-doped crystals 
demonstrating that the key to the CMR effect is to be 
found in the non-metallic regime. The ionic size of the di- 
valent substituent exerts its effect on magnetic and elec- 
tronic properties through local tilting of the oxygen octa- 
hedra [Q and the concurrent bending of the active Mn-0- 
Mn bonds. This inhibits the formation of metallic bonds 
and leads to charge localization, polaron formation, and 
possible charge segregation [||j9j. Evidence in favor of a 
percolation picture comes from a Monte-Carlo simulation 



of a random field Ising model that assigns conductivity 
zero and unity to bonds between neighboring antipar- 
allel and parallel spins respectively [Q. It both pro- 
duces CMR and emphasizes the importance of random- 
ness. Experimental evidence was provided by Jaime et al. 
[ pT| who extracted the field and temperature dependent 
metallic-bond concentration c(if, T) from the resistivity 
via the effective medium approximation, and showed that 
it also describes the thermoelectric power data. Direct 
evidence for coexisting polaronic/insulating and metallic 
components have been reported from neutron [ p^ , [ p^ , 
electron 14 1 and Raman scattering ||l^ studies. 

Discussions of the CMR effect in percolation terms 
]l6|-p^ have generally neglected a central point: that the 
percolating entities are formed thermodynamically as the 
temperature is lowered, nucleated by the intrinsic ran- 
domness of a doped material and amplified by the ten- 
dency for polaron formation and charge segregation . 
Griffiths first pointed out, in the context of dilution 
(bonds randomly assigned values of or J), that singu- 
larities arise in thermodynamic properties in the temper- 
ature range TjJ?"'^ < T < Tq between the random tran- 
sition T^^"'^ and the "pure" transition temperature Tq. 
Bray |2C|j2l[| extended the argument to any bond distri- 
bution that reduces the transition temperature, terming 
the regime between T^'^"'^ and Tq the "Griffiths phase." 
In this paper, we use the Bray model to demonstrate 
that the zero-field transition is an example of a Grif- 
fiths singularity. To treat the field dependence, we con- 
sider two possible exchange energies: Jmeti associated 
with the double-exchange process on metallic Mn-O-Mn 
bonds, and Jms, the residual ferromagnetic interaction 
that exists on insulating bonds. The proportion of metal- 
lic bonds changes with field and temperature, monitored 
through the resistivity-based metallic fraction c{H,T). 

Among the effects that characterize the Griffiths phase 
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are a distribution of susceptibilities and slow spin dy- 
namics. The latter have been observed by muon re- 
laxationj 22p3] and show up strongly in noise measure- 
ments |24]. We focus in this work on the nature of the 
susceptibility and magnetization, and on the behavior of 
the heat capacity in applied magnetic fields, tying the 
last to the CMR effect itself. Figure 1 shows the in- 
verse susceptibility x~^(T) — H/M{T) of a single crys- 
tal sample of Lao.rCao.sMnOa, where M{T) is the mag- 
netization in a field of 1 kOe. The high-temperature 
behavior is of Curie-Weiss form, with a paramagnetic 
Curie temperature 9 = 202 K. Well above Tc = 218 
K the slope of the experimental curve corresponds to 
5 w 8, much shallower than the Curie-Weiss law for 
spin Seff — 1.85 (the weighted average of 5 = 3/2 and 
S' = 2 appropriate for this sample), shown as a dashed 
line. Further, there is a sharp downturn in x^^iT) before 
the Curie temperature is reached. This alone identifies 
the transition as a Griffiths singularity, characterized by 
a susceptibility exponent less than unity ||2^; that is, 
X~^{T) oc (T — T^^"^)^"'*', with a non-universal positive 
exponent A < 1. 

In his treatment of the Griffiths phase. Bray [ p0|j2l| 
argues that the distribution of cluster sizes can be char- 
acterized by the longitudinal susceptibility matrix Xij — 
g^f^B i^i^j) I^bT or rather by its inverse, whose eigen- 
values ii/C are distributed according to where 
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Here C is the Curie constant of the material and each 
H represents an effective Curie- Weiss temperature dif- 
ference T - SeffiT). The Qeff{T) may be equally 
well related ||2^ to the eigenvalue distribution of the 
exchange matrix which the double-exchange mechanism 
makes temperature and field dependent. The close con- 
nection between the effective exchange energy and the 
metallicity of the corresponding Mn-O-Mn bond ties the 
Griffiths behavior to the CMR effect. Bray ignored the 
negative-power prefactor, but it is required if is to 
decrease for large fj, and to rise sharply near fj, — as 
jrand approaclicd |2^. As all clusters are finite, none 
is fully ferromagnetic so we must have p{fi — > 0) = 0. 
However, there is a pile-up of large clusters as the tran- 
sition is approached, leading Bray to deduce that A{T) 
vanishes as, 



A{T) = To 
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where Tq is a parameter and [3 = 0.38 is the order- 
parameter exponent for the pure system, assumed to be 
?>D Heisenberg-like. Because fi = T for free spins, we 
cut off the distribution and calculate the average 
inverse susceptibility (/i) / C, using 



The resul is shown as a solid line in Fig. 1, where the 
Curie constant is 

N9^f^lS,ff{Seff + l) 



C 



and TJ'''^'^ 
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0.074 K 
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(J — 223 K. The free parameters are the charac- 
teristic temperature Tg = 9.5 K and the non-universal 
prefactor exponent, 7 = 1.59 ±0.5. The calculated curve 
is very close to the power law x~^{T) cx (T - r^and)0.72^ 
so that A = 0.28, as seen in the inset to Fig. 1. The 
values of Tq and 7 are strongly correlated; the data can 
be fit with values of 7 in the range cited, resulting in 
values of A in the range 0.2 < A < 0.3. Analysis of the 
quantum critical point of UCu4Pd in terms of a Griffiths 
singularity yields the much larger value A = 2/3 ||2^ . 

The downturn in the x~^{T) data close to the tran- 
sition temperature is more abrupt than Bray's approach 
predicts and J^J?""^ is somewhat higher than Tq =218.2 K, 
the location of the zero-field heat capacity peak shown in 
Fig. 2 1^ that we take to mark the true transition tem- 
perature. Griffiths, in his original paper suggested 
that the susceptibility would tend to diverge in advance 
of the onset of true long range order, and that may ex- 
plain the present result. The nonanalytic behavior of 
the magnetization at small magnetic fields is readily ob- 
servable in the vicinity of Tc where the magnetization 
rises so abruptly that, considered in terms of conven- 
tional critical behavior M oc H^^^ , it requires S « 13, far 
from any standard value. Magnetization data taken at 
218 K, showing this rapid rise, are shown in Fig. 3. The 
upward curvature in the x~^i'^) data at high tempera- 
tures has also been reported by, for example, de Teresa 
et al. [p6| , who find a linear regime at high temperatures 
with a slope consistent with the Curie constant in Eq. (|^) 
and an effective Curie- Weiss temperature w 1.25Tc; it is 
tempting to assign this to Tq- 

Figure 2 shows the heat capacity of the sample in 
fields of T, 3 T and 7 T after subtraction of a smooth 
background that fits the data at temperatures well away 
from the peak. The sharpness of the zero field data 
might signal a first-order transition, but there is no evi- 
dence of hysteresis. The shift in the heat capacity peak 
is uncharacteristic of ferromagnetic transitions, as dis- 
cussed previously To treat this behavior, we em- 

ta modified version of the Oguchi cluster method 
a simple improvement to mean-field theory. In this 
method, the energy of a pair is treated exactly, while 
coupling to its z — 1 neighbors is handled in mean-field 
theory. The magnetic part of the double-exchange en- 
ergy is xt[cos{9/2) — cos(0/2)], where 9 is the classical 
angle between core spins, t is the transfer energy, x is the 
doping level and the bar denotes the angular average. 
For the Oguchi calculation, we use the quantum analog 
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where 5*0, the total spin of the two manganese S'c = 3/2 
cores and the shared eg electron, is in the range 1/2 < 
Sq < 7/2. The pair is coupled to the magnetization of its 
z — 1 neighbors via either metallic Jmet or insulating Jins 
exchange bonds. Assuming spherical clusters (of undeter- 
mined size), we follow the effective-medium approach of 
Jaime, et al. and extract the metallic bond concentra- 
tion c{H, T) from the measured resistivity p{H, T) (Fig. 
2a) by extrapolating the fitted the high-temperature (ac- 
tivated) and low temperature (power-law) resistivities. 
The resulting c(iJ, T) is shown for several applied fields 
in Fig. 2b. The effective field acting on the pair is, then, 

Heff = H + 2{z - l)SeffJeffm{H,T)/gfiB, (6) 

where m(H,T) is the reduced magnetization, to be cal- 
culated self-consistently, and Jeff — c{H,T)Jmet + [1 ~ 
c{H, T)]Jins- The resultant m{H, T) values are then used 
to determine the pair energy (EdeiSo)) within the Oguchi 
approach. 

The field and temperature dependent heat capacity 
calculated from {Ede{So)) with a single amplitude pa- 
rameter, shown as solid lines in Fig. 2c, is in good 
agreement with the data. The insulating exchange en- 
ergy is determined from the observed Curie constant to 
be J^ns = 3kBe/2zSeff{Seff + 1) = 0.82 meV. We ob- 
tain the transfer energy t from the spin-wave stiffness D 
according to t « (25*^ -I- l)D/xa'^ |29). The energy differ- 
ence between parallel and antiparallel configurations is 
Ede{7/2) — Ede{i/2) = 3xi/5; for a Heisenberg interac- 
tion —2Jmet~§i-~§2 the difference would be 4,3^ Jmet giv- 
ing Jmet = 3x^/205'^. For D = 150 mev A^, which is com- 
mon to most manganites at this doping level we find 
t = 130 meV and Jmet = 1-46 meV. Clearly, the back- 
ground subtracted from the data has ignored magnetic 
contributions that remain well below Tc- The shoulder 
in the zero-field curve arises from the tail in c(0, T) in the 
vicinity of Tq- This may signal a break-down in the ef- 
fective medium approximation, as very large clusters are 
known from noise measurements p^j30| to have a dra- 
matic effect on the zero-field resisitivity near the transi- 
tion. The self-consistent values ■m{H, 18 K), shown as 
open circles in the inset of Fig. 3, match the measured 
critical isotherm. Neither Jmet nor J^s determines the 
location of the zero-field heat capacity peak at 218.2 K. 
Rather, the rapid buildup in both the size and number 
of magnetic clusters that give rise to the Griffiths behav- 
ior of the susceptibility causes a concomitant increase 
in the conductivity of the material. As in Griffiths' di- 
lution picture, the existence of Jmet serves to pull the 
transition temperature upward from the values expected 
from Jins- Insofar as it is characterized by the spin- wave 
dispersion at low temperatures, Jmet is essentially the 
same for all CMR materials. Therefore, the stronger 
tendency toward self-trapping that occurs as the Mn-0- 
Mn bond angle is reduced, lowers Jins relative to Jmet 
and increases the dominance of the Griffiths phase. 

We have argued that the CMR effect should be treated 



in the context of a Griffiths singularity driven by intrinsic 
randomness, the combined effect of doping, the tendency 
for charge segregation, and the self-trapping effect asso- 
ciated with polaron formation. Bray's treatment of the 
distribution of susceptibilities explains the unusual down- 
turn in the measured susceptibility and the magnitude of 
the effective paramagnetic moment. More importantly, 
we have used the Griffiths idea of a distribution of ex- 
change energies to demonstrate that the magnetic tran- 
sition is driven by the accumulation of larger exchange 
energies on metallic bonds as the system begins to order. 
Magnetic fields align large clusters, turning bonds along 
their shared boundaries from insulating to metallic. This 
has the effect of shifting the effective percolation point 
upward in temperature, and connecting the conductive 
network at temperatures well above the zero-field per- 
colation point. This is the essence of the CMR effect. 
A modified Oguchi model that treats pairs in the con- 
text of the double-exchange mechanism, and that uses 
the resistivity-derived, metallic-bond concentration, ac- 
curately captures the nature of short range correlations, 
although not the effect of larger clusters. Future work 
will address the self-consistent calculation of the concen- 
tration and magnetization without the need for empirical 
input from the CMR data, and will apply this analysis 
to other CMR transitions. 

We gratefully acknowledge the contribution of Y. 
Tomioka and Y. Tokura in providing the sample used 
for this study and have benefitted from discussions with 
M. B. Weissman. This work was supported in part by the 
U.S. Department of Energy through contract DEFG02- 
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FIG. 1. Ratio of the magnetic field H to the magnetization 
M at 1 kOe. The dashed line is the expected Curie- Weiss be- 
havior in a 1 kOe field for spin S = 1.85 and Tc""'^ = 223 
K while the solid line is the result of the Griffiths-phase, 
H = calculation. The inset shows the Griffiths phase re- 
sult vs reduced temperature on a log- log plot, demonstrating 
its power-law behavior. 



3 



300 



-P 150 



o 100 
CO 
Q. 
CD 

^ 50 
CO 

o 
X 



a) 




. T 




■ 3 T 




' 71 







100 ISO 200 2S0 

T(K) 




100 150 200 250 300 350 

Temperature (K) 

FIG. 2. Heat capacity of a Lao.rCao.aMnOs single crystal 
at several fields. The solid lines are computed from the mea- 
sured metallic-bond fraction using the Oguchi model. The 
resistivity and metallic fractions at the same fields are shown 
in insets a) and b) respectively. 
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FIG. 3. Magnetization vs applied field at 218 K. The inset 

is a linear plot; the main panel is logarithmic, showing the 
small slope given by l/J. The open circles in the inset are 
calculated at 218 K within the Oguchi model. 
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